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Abstract: In this work we introduce a new polynomial representation of the Bernoulli numbers 
in terms of polynomial sums allowing on the one hand a more detailed understanding of their 
mathematical structure and on the other hand provides a computation of B 2n as a function of 
B 2n _2 only. Furthermore, we show that a direct computation of the Riemann zeta-function and 
their derivatives at k e Z is possible in terms of these polynomial representation. As an explicit 
example, our polynomial Bernoulli number representation is applied to fast approximate compu¬ 
tations of C(3), C(5) and ((1). 
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1 Introduction 

Much was written on Bernoulli numbers |ffi[2,[3, 4^ [5] since Jacob Bernoulli has discovered these 
fascinating quantities f6j. Their importance is widely documented as they appear in a variety 
of scientific fields of importance EE, 8], [9J . In particular, their close relationship with the Rie¬ 
mann zeta-function [ITOl ITTl iT2ll . shows the need for a better understanding of their mathematical 
structure. Bernoulli numbers are rational, their fractional part is known by the results of Karl von 
Staudt and Thomas Clausen, often cited as the Staudt-Clausen theorem Ifl3l fl41. Furthermore, 
several direct representations exist EQj, El El, all show up as complicated expressions in form 
of double-sums providing not much information on the enumerator as well as on the denominator 
of these numbers. No simple rules to compute Bernoulli numbers had been established so far and 
in consequence explicit formulas are more or less of pure academic interest [E| • On the other hand 
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simpler representations could be of great practical importance, for example, in the field of analyt¬ 
ical number theory. For a direct computation of the Riemann zeta-function and their derivatives 
for all k G Z in terms of Bernoulli numbers is possible. Let us remark first on these relations. 
Some of such relations are history know like the famous Euler representation 12311241 

a2n) = , <U) 

or for negative integers Ifl 2 l 


C(-n) 


n + 1 


( 1 . 2 ) 


More delicate is the direct computation of (’(2n + 1) and of ( 7 (2n). Explicit relations were found 
by us and published in 1994 IfTTTl and similar one later by Adamchik et al. Ifl2ll25l . Some remarks 
on the corresponding computational scheme will be given with an explicit example later on, where 
we demonstrate the direct application of our polynomial Bernoulli number representation to ( (3). 
At first, we found 


and 


c(2n+i) = (—r + %n 


2(2tt 


\2 n 


(2 n) 


C'(2n) = (-)’ 


( 2 yr) 2n „ / 2 n r . . 

:-£>2n ( —— L2n-1 ~ 7 E ~ ln{2’K) 


2(2 n)\ 


B- 


2 n 


(1.3) 


(1.4) 


where 7 e denotes the Euler constant. The L-numbers were first introduced by Bendersky in 
context with the logarithmic Gamma-function Il26il27ll . The first derivative on the odd numbers 
C'(2 n + 1 ) may be computed as follows 


C'(3) 


2tt 2 C"(-2) + 



2 + C(3) , 


(1.5) 


where 


c"(— 2 ) = 3 L 2 


17 B2n+2H2n-2 

108 + (2n — l)( 2 n)( 2 n + l)( 2 n + 2 ) 


(1.6) 


can be written in terms of a Dirichlet series with the Bernoulli numbers and the harmonic numbers 
If28l involved. Furthermore, we recall one of our earlier results for ( (3) published in IfTTI where 
we have introduced a fast converging series representation based on Benderskys L-numbers ll26l 


C(3) 




C(2n) 

2n(2n + l)( 2 n + 2 ) 



(1.7) 


Using the well-known Taylor-McLaurin series representation for C(s) ll29l 
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with sG — N, one finds for the first derivative at s=-2 


C'(-2) 


1 

36 



B 2 n 

3)(2n — 2)(2n 


l)2n ' 


(1.9) 


Following Bendersky If26l the logarithmic gamma function lnr n (x) of order n may be written as 
follows 



A n (x + 1) = x n ln(x ) — 

L n + ItiY' ni^c) 

( 1 . 10 ) 

where, for example, A 2 ( 

x + 1 ) is defined as ll26l 



A 2 (x + 1) 

(l 3 l 2 l A 

= -X + -X + -x 

V3 2 6 J 

ln{x ) 



00 

- r 3+ ^ 2+2 E 

n=2 

Q 2n+l 

-D 2n +2 ^ 

( 1 . 11 ) 


(2n — l)(2n)(2n + l)(2n + 2) ' 

At x = 1 it follows 

A 2 ( 2 ) = C'(-2) 

= l 2 . 

( 1 . 12 ) 


The L-numbers had been computed by Bendersky for all integer numbers n If26ll . For n=2 we 
have 


L 2 = — I - -inf- 
48 l 2 V 3 




C(2 n) 


(2n)(2n + l)(2n + 2) V6 


n =1 


2 n' 


(1.13) 


This procedure allows a direct computation also of ('(—n) as shown by us before ifTTI . 

C(~ n ) = ~L n + 


B n +1 


y-- 

n + 1 q 
q =1 


(1.14) 


Within this procedure it is possible to find fast converging series in terms of the L-numbers for all 
C(2n-1) values. For example, for n=3 and n=5 one gets 


C(5) 



n= 1 


257T 4 7T 4 / 7T \ 

-1- In ( — ) 

12528 1044 V 3/ 

_CM_ 

2n{2n + l)(2n + 2)(2n + 3)(2n + 4) 



(1.15) 


and 


C(7) 


72tt 2 

"659" 


C(5) 


32tt 6 

5931 


£ 

n= 1 


2tt 4 497t 6 7r 6 /7r\ 

"~ 197y ^ + 5337900 + 266895 ^ V37 

_CM_ 

2n(2n + l)(2ra + 2)(2n + 3)(2n + 4) (2 n + 5)(2n + 6) 



(1.16) 


3 
















For related constants like U 2 or U 4 we found: 


and 


\n+l 


U2 ^ ^ (2n — l) 2 

OO 

*E 


n= 1 

= 1 - 


n+l nL 2n (K\ 2n 


n= 1 


(2m)! V2 


u* = E< 

n =1 

9 

2 ” 


\n+l 


(2n-l) 4 


7 7 4 00 

2 /n(2) - -U 2 - -C(3) + -£( 


(1.17) 


,+ l ^ 3 ^ 2 n +2 ( 7 T \ 2 "+ 2 


n=l 


(2n + 2)! \2/ 


(1.18) 


These examples show that not only zeta values can be expressed by Bendersky’s L-numbers and 
as a consequence by the even Bernoulli numbers. With this work we introduce a new polyno¬ 
mial representation of the even Bernoulli numbers which allows for a fast computation of these 
constants, which is comparable with the computation by use of available BBP-like formulas. Fur¬ 
thermore, our method allows for the use of polylogarithmic identities on which BBP formulas are 
typically based. 

The paper is organized as follows. In the next section we introduce a new formula for the even 
Bernoulli numbers which serves as a basis for our polynomial representation. We present a prove 
of this formula discuss it in context with other well-known series representations of Bernoulli 
numbers, giving a special emphasize on the numerical effort which is needed to compute B 2n . In 
section 4 and 5 we present first applications of our formalism to the computation of £(3), C(5) 
and £(7) and compare with corresponding BBP-type formulas concerning the numerical effort. 
The final section 6 gives a summary of the results and a short outlook. In the appendix section 
we present as a further application of our polynomial representation the computation of B 2 „ as a 
function of B 2n _ 2 only. Based on this formula an even faster computation of zeta numbers by one 
order of magnitude is possible. 


2 An alternative formula for the even Bernoulli numbers 

In this section we present the basic equation from which the Bernoulli numbers B 2n can be com¬ 
puted. 

Theorem 2.1 


It follows for B 


2 n 


Bo n — 


2 2n+1 (2 n)\ 
2 2n - 2 


E 


\fc+i 


n + l 


—' (2n + k)\ y k + 1 


Eh' 


1=0 


k 

l 


- ir 2n 
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( 2 . 1 ) 
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ma 
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This work 

Bio 

2.5 

1.5 

1.5 

3.0 

2.0 

1.0 

B 20 

5.0 

1.8 

2.0 

5.2 

3.0 

1.0 

B 30 

4.8 

1.6 

1.7 

5.3 

3.1 

1.0 

B 40 

4.8 

1.6 

1.8 

5.4 

3.1 

1.0 


Table 1: Runtime comparison between different recurrence relations for Bernoulli numbers with 
our present formula. Shown is the factor which results from an explicit calculation of the even 
Bernoulli numbers for n=10, 20, 30 and 40. 


The formula shown above represents an alternative computational scheme for the Bernoulli num¬ 
bers, as only the even numbers are calculated, where the outer summation index runs up to n 
only. All other formulas available from literature E [T5j [T6. 17] including the representation 
which was introduced in 2009 lf30l compute the Bernoulli numbers B n by a summation up to n. 
A direct computation of the Bernoulli numbers by use of Eq. (2.1) is possible, where a runtime 
comparison with a variety of other existing recurrence formula reveals that our formula allows 
for the computation of B 2n numbers which is typically faster by factors ranging between 1.5 
and 5.4. Corresponding numerical results are presented in Table I: Furthermore, and this is the 
most important point, our double-sum like representation allows a direct link to the so called 
central factorial numbers described in detail by Riordan If22l . which finally offers the possibility 
to introduce a polynomial representation of the Bernoulli numbers not known before. We show 
later on that this polynomial representation allows for a fast computation of all odd zeta-numbers 
£(2 n + 1) and related constants discussed in the introduction. 

To prove the theorem we first list some helpful identities. 


k — 21 I i + 1 \ ( k — l 
i + 1 1 k — l 


k - l - 1 


k — l — l 
l 


( 2 . 2 ) 


This identity is simply proved by a direct calculation of both sides. 


71—1 

E 

i =0 


k-2l 


n + 1 — 


k - 21 


n 

i 


n 


Ti 4“ 1 — k 4- 21 


(2.3) 


Identity (2.3) is found in Il22l I3TT1 . 
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77+1 


n 


7i + l — k + 21 \ k — l 


k-l 

l 


77+1 

k-2l 


n — k + 21 

l 


(2.4) 


Again, identity (2.4) is proved by a direct calculation of both sides. Now we start our proof 
with the well-known recursive formula [O]: 


B 2 n — ~ — 


277 + 1 


n— 1 

E 

1=0 


277 + 1 
21 


B - 


21 ■ 


(2.5) 


This recursive relation can be used to formulate an iterative solution for the Bernoulli numbers 

n— 1 


Si+i( n ) = 

k=i 

i=l, 2, 3, ...G N where it follows for Si( 77 ) 

n —1 

s i( n ) = 


2n + i + 1 
2k + i 


Si(k) 


( 2 . 6 ) 


1^-4/ 277 + 1 


7=1 


21 


= -( 2 2n — 277) 


(2.7) 


In second order this recursive formula gives 

S 2 (n) = ^3 2n - (277 + l)2 2n + 2 t7 2 + 77 - | 

O O 

Then the Bernoulli numbers B 2 „ result in 

2 B 2n 1 2^(77) , 25 2 (77) 


( 2 . 8 ) 


277 


+ 


2n 277(277 + 1) 2ti(277 + 1) (277 + 2) 

2S 3 (n) ± 


(2.9) 


277(277 + 1) (277 + 2) (277 + 3) 

For a direct computation of (2.9) we convert in a first step all partial sums Sj(t 7), i=l, 2, 3,... 6 N 
to a more appropriate form. For 7 = 1,2 and 3 it follows explicitly 


Si( 77 ) = ^2 2n+1 - i(277)l 2n+1 


S 2 (n) = l3 2n+2 +-2 2n+2 - -l 2n+2 - 2(277 + 2)^(77) 

8 4 8 


( 2 . 10 ) 


-(2t7 + 1)(2t7 + 2) , 


^ 3 ( 77 ) = J_4 2n + 3 + I 3 2«+3 _ I 2 2n+3 _ 3 1 2n+3 _ 3 ( 2n + 3 ) 5 2 (t 7 ) 
16 8 8 8 

- 3(277 + 2)(2n + 3)5i(77) - ^(2n + 1)(2ti + 2)(2ti + 3) . 


( 2 . 11 ) 


( 2 . 12 ) 
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Formula (2.12) represents a ternary sum because of the repeated use of the recursion formula 
(2.6), and as a consequence S n {n) denotes a n-fold sum. At a first glance, this fact sounds very 
discouraging but fortunately most of the terms cancel if all partial sums Sfn) with index i=l,...,n- 
1 will be inserted in the corresponding n-fold partial sum S n (n) with index n. To illustrate this 
procedure we compute the partial sum with index i=3 

Example 2.1 


iiis 

C _ __j2n+3 I _o2n+3 _^2n+3 _i2n+3 

31 J “ 16 8 8 8 

- 3(2n + 3) (- 3 2n+2 + ~ 2 2n+2 - il 2n+2 
+ 3(2n + 2)(2n + 3) Q 2 2n+1 - i(2n)l 2n+1 

- i(2n + l)(2n + 2)(2n + 3) . 


(2.13) 


This algorithm has to be applied to each partial sum Sfn) in a recursive sense. The result for 

Si(n) is 


Lemma 2.1 


with 


s>(n) = (fc-oi 

k =0 1=0 


i + 1 — k 


k—l 


2 n + i 
k-l 


i + 1 — k 


(2.14) 


MO = 


7 Tl Til 

cos - sm — 

2 2 


7 Tl Til 

cos -h sm — 

2 2 


i + l — k 
ti{l) 

i + l — k 
t 2 {l) 


(2.15) 


with ti(l) = | j + ^-) and t 2 (l) — f — | — The index t i(l) produces the numbers 0, 0, 
1, 1,2, 2, 3, 3, for 1=0, 1, 2, 3,..., k. Thus has the same values for 1 even or odd. The second index 
t 2 (() is slightly different as this sequences starts with -1 for 1=0. The proof can done by induction. 

As an example, we show in table II the complete list of numbers which result for S l (n), where 
i=l,...,6 and n=l,...,10: 

Proof. 
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Si(n) 

S 2 (n) 

S 3 ( n ) 

S 4 (n) 

S 5 (n) 

S 6 (n) 

n=l 

1 

0 

0 

0 

0 

0 

n=2 

6 

20 

0 

0 

0 

0 

n=3 

29 

392 

1680 

0 

0 

0 

n=4 

124 

5112 

73920 

369600 

0 

0 

n=5 

507 

55220 

2000856 

30270240 

16816800 

0 

n=6 

2024 

544700 

43099056 

1462581120 

2.287+ 10 

1.372+ 11 

n=7 

8185 

5135184 

821292576 

5.475+ 10 

1.788+ 12 

2.868+ 13 

n=8 

32760 

47313584 

1.459+ 10 

1.772+ 12 

1.062+ 14 

3.384+ 15 

n=9 

131063 

430867484 

2.489+ 11 

5.250+ 13 

5.365+ 15 

2.991+ 17 

n=10 

524278 

3900612564 

4.140+ 12 

1.470+ 15 

2.444+ 17 

2.218+ 19 


Table 2: Computed numbers for Si(n ) with i running from 1 to 6 and n running from 1 to 10. 


Using now (2.9), doing few elementary manipulations, one ends up with a first expression for 
B 2n in form of a ternary sum only 


n i ti(k) 

B 2n = 2 2n {2n)\Y J Y,H^ k+l+l 

i=0 k =0 1=0 


k — 21 


i — k + 21 

l 


2n + * — 2k + 4/ + 1 
i — k + l + 1 

/ 2-)_i_^ \ 2n+z—&+2Z+1 

V 2 ) _ 

(2 n + i — k + 2/ + 1)! 


(2.16) 


Making use of lemma (2.1) the different terms in (2.16) can be rearranged, leading to a simpler 
formula for B 2n 




2 2ri 


i=l fc=0 i=0 

« \ I (^) 


2_|_][_ fc \ 2t 7/—|—z — fc+2/ + l 


k — l + 1 
A further rearrangement results in 


(2 n + z — k + 21 + 1)! 


— 


o2n/n„ \ | n 


2 2n- 


n 

k-l 


k =0 t=0 


k-l\ (Sfi) 


\ 3n—fc+2/ 


l / (3n-fc + 2Z)! 


71—1 

E 

i=0 


k - 21 


-l 


n 


(2.17) 


, (2.18) 


where we have interchanged the sum over i with the other two sums. As the variable i appears in 
the inner sum only this offers the possibility to reduce the ternary sum to a binary one. This is 
done by use of lemma(2.2) lf22l[3TTl . 
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k — 21 

As the contribution from the Binomial coefficient | I is zero it remains for B 2n 


n 


r)2n( n \ | n-Hi(k) 


n + 1 
k — 2l 


k =0 i=0 

n — k + 21 A 


_£> \ 3n—fc+2Z 


/ I (3n-k + 2l)\ 
where we made use of lemma (2.3). Next, we rewrite (2.19) as follows 

,„ + i 2 2 ”(2n)! 


= (-) 


and use the relation 


E 


2 2n ~i - 1 ^ (3 n-k)\ 

k =0 v ’ 


t\(n—k) 

* E (-)' 

1=0 


n — k 
l 


n — k 


n + 1 
k 

3 n—k 


(2.19) 


( 2 . 20 ) 


tl(n-fc) / 

E (-)' 

i=o \ 


n — k 
l 


(n — k — 21) 


3 n—k 


ti(fc) / 

EH- 
1=0 \ 


(k-2lf n+k , (2.21) 


l 


which is simply proved in replacing n-k by k as the summation is symmetric in these indices. 
Finally we result in 


B 2n — 


2 2n (2 n)\ 


2 2r 

ti(k) 


rzrE( 


\fc+l 


n+ \ \ (|) 


1 \ 2n+k 


k=0 


S-' (;) 


n — k J (2 n + k)\ 

(k - 2lf n+k 


( 2 . 22 ) 


Finally the summation over ti(k) can be replaced by a summation over [|] because the inner 
sum behaves symmetric with respect to the index ti(k). This last step completes the proof of the 
theorem as the binomial coefficients in (2.1) and (2.22) are identical. □ 


3 The a-polynomial as a generator to compute Bernoulli num¬ 
bers 

We define the inner sum in (2.1) as a polynomial of order n in the variable k e N with coeffi¬ 
cients a\ n \ 1=1,...,n, which can be directly computed in a numerical sense by use of the definition 
shown below. To define the inner sum as a polynomial at a first glance seems to be quite uncon¬ 
ventional, but this definition is based on corresponding relations which can be found, for example, 
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in lf22ll3TI for k=0,l and 2. The generalization to k>3 is straightforward, where explicit formulas 

(n) 

for the coefficients a\ will be given at the end of this section. We have then 

Proposition 3.1 


n 



1=1 


k l 


2(2n)! y ^(_)j ( ^ 

( 2 n + k)\ yj 


,k 

'2 


j) 


k-\-2n 


(3.1) 


Therefore, it follows for 5 2n 


B2n ~ 


2 2n 

2 2n - 2 


Eh ‘ +1 

k=1 


n + 1 
k + 1 


E 

z=i 


(3.2) 


Next, we show that the o-polynomial can be defined for negative integer numbers. Furthermore, 
it will be demonstrated that a computation of the o-polynomial for negative integers provides 
a new option to define the Bernoulli numbers. We start with a recursive formula which the a- 
polynomial fulfills. It is 


Lemma 3.1 


i=i 


k(k - 1) 

(2 n + k){2n + k — 1) 


E of (*-*>’ 

i=i 


1 j 2 2n(2n — 1) ( n -i)u 

4 (2n + k)(2n + k — 1) ^ 1 


(3.3) 


with Oq = 1. 

To prove the above lemma let us first recall the following recursive formula obtained by Rior- 
dan El, for the so called central factorial numbers 0 


T(n,k ) = T(n — 2,k — 2) + -k 2 T(n — 2,k), (3.4) 

with 

^)^Eh'(:)(H"- (3 - 5) 

From this it follows immediately 


T{2n + k, k ) 


2 n + k 
k 


E“l n) *' 


i=i 


(3.6) 


1 Riordan discusses these numbers in context with the so called central difference operator 5 in chapter V of his 

book ’Combinatorial identities’. 
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Now it is easy to complete the proof. First we write 



(3.7) 


Substituting here (3.6) the recursive formula for the o-polynomial follows directly for ke N. 
Furthermore, we have shown that the central factorial numbers are closely related to the a- 
polynomial, which can be regarded as a key quantity in computing the Bernoulli numbers. 

Using the notation 

Definition 3.1 


n 



(3.8) 


i=i 


it follows for example for k=l 


Example 3.1 



(3.9) 


This way the sum /l (n) (k) may be calculated for all k e N. To extend the calculational scheme 
to negative integers we introduce the following relation 


Lemma 3.2 



(3.10) 


To prove the above relation we use the following identity ll22l[31 'j 


Proposition 3.2 



(3.11) 
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Proof. 


Writing now 


tc&H-kr = 

m= 1 \ / 1=1 m= 1 

= fcl n + k n IT 


k 


m= 1 Z=1 


k + l \ l 


E “£>(-*0” ■ 


m=l 


(3.12) 


□ 


As a consequence we find the result that A^ (—1) is ’’proportional” to B 


2 n 


Example 3.2 


o2 n o 

= - - . 

Using (3.3) it follows further 

j4 W(_ 2 ) = (2n — 1) B 2n , 


and for example 

i o2n _ o 

= -(2n)(2n-l)-^B 2 „ 

1 2 2 "- 2 - 2 
- i (2n-l)(2n-2)^ z ^5 2n _ 2 , 


(3.13) 


(3.14) 


(3.15) 


where in the expression for A^ n ^(—3) both B 2n and B 2n -2 appear. More detailed information 
about the properties of the A (n) (k) can be obtained by changing to an alternative representation, 
which was first introduced by Riordan lf22l . again in context with the so called central factorial 
numbers 


Lemma 3.3 


A {n \k) 


1 

2^ 


2 n + k 
2 n 


with a n , 0 = 5 n ,o and a n ,i=l . 


n 


^ ^ dn,l 
1=0 


2 n + k 
2 n + l 


Proof. 


(3.16) 
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Following Riordan ll22ll we first have 


2 2n T(2n + k, k ) 


n 


z=o 


2 n + k 
2 n + l 


Using (3.6) it follows 


T{2n + k , k) 


1 

2^ 



2n + k 
2 n + l 


2n+ M 4<”>(fc) 

2n I 


(3.17) 


(3.18) 


and this completes the proof. □ 

The coefficients a n j are defined by the recursive formula presented below Il22ll 


^n,l 


n— 1 

^ ^ 1 

z=Z—1 


2n + / — 1 
2n — 2i 


(3.19) 


with a„ i=ao,o = 1 and a n 0 = 0 for n>0. It should be mentioned here that no obvious connection 
between the central factorial numbers and the Bernoulli numbers exists. This is because the re¬ 
cursive relation (3.19), which allows to compute the central factorial numbers is linked by (3.16) 
only to the a-polynomials but not directly to the Bernoulli numbers. Thus, it becomes clear that 
the Bernoulli numbers were not discussed by Riordan in context with the central factorial num¬ 
bers. As an example it follows for n=2 


Example 3.3 


A^ik) 


a 


^k + offik 2 


1 

¥ 


k + 4 
k 


i k2 + \ Bik • 




k + 4 
6 


(3.20) 


where the Bernoulli number B 4 appears as the coefficient of the lowest order polynomial term. 
Computing (3.16) for k = -1 we find for B 2n 

1 n 

B ’ 2 n = 22 n _ 2 ) + an ’ 1 

1=1 

In a next step we substitute for a n j 



(3.21) 


Definition 3.2 
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<v< = ( 2n 2 t ‘) < 3 - 22 > 

by introducing a new type of polynomials. The result for B 2n is then 

3* = ( 2 g» 2 -2)6“ E (-)' +11 P< " +1_ ' , (") ■ 0.23) 

The P-polynomials can be obtained from (3.22) but they are also connected to the ^-polynomials 

Lemma 3.4 


p(n+ 1 -0( /l ) = 


2ncn( \l 


2 2n 6 


/(2n)! 


E 

fe=i 




and vice versa for A^ n \k) 


f M pOH->-0(„) . 


2 2 n 6 n 

Z=1 


l J 


(3.24) 


(3.25) 


Proof. 

First, the relation for A^fk) results by substituting (3.22) in (3.16) and by reformulating 
the product of the three binomial coefficients. The above formula for the P-polynomials is then 
proved by substituting (3.24) in (3.23). This gives 

o2n n [ ( 7 \ 

Bln = E E(-)‘ +1 k ' < 3 - 2 ® 

1=1 k =i V / 

Replacing now l by n + l — 1 it follows 

o2n n n f 1 ~ l / r, _l_ 1 _ / \ 

Bin = E • < 3 - 27 > 

z=i fc=i V / 

The double sum can be written as 

E (-) t+1 ( l ) AM(k) + ... + E (-) t+1 ( l ) A(n> (k ) 

= D-)‘ +1 f H J) . < 3 - 28 ) 
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which finally gives (3.2). □ 

As an example, we have for the first three P-polynomials 

Example 3.4 


p(i) (n) = I pW(„) = A, ^ 7 43) . (3.29) 

All P-polynomials can be regarded as ordinary polynomials, apart from the first one. 

It was shown by (3.10) that the o-polynomials can be defined for negative integers. The cor¬ 
responding formula computed as a function of the P-polynomials is given by 

Lemma 3.5 


A^ n \-k) = (-)’ 


(2 n)\ 
6 n 2 2n 




\i +1 


«=i 


n + A: — l 
k 


p(0( 


n 


(3.30) 


Proof. 


Computing (3.25) for k = -1 gives 

(2 n)\ 


A^ n \-k) = 


6 n 2 2r 


Eh ' 1 


z=i 


l + k — 1 


P (n+1 " z) (n) . 


Replacing again l by n + l — 1 it follows 

= MD-r'An + lV ” + *"' 


t=l 


n + 1 — / 


p(0( 


n 


6 " 2 2 


i=i 


fc - 1 


(3.31) 


(2 " )! £(-)'«(„ + fc - /) f " + ( " \ ~ 1 ] P<‘>(n) ■ 


(3.32) 


□ 

The relation (3.32) allows an explicit definition of the a \ n ' > -coefficients by comparing directly 

(n) 

the corresponding coefficients in the polynomials on both sides. For a\ we found the interesting 
result 


Corollary 3.1 
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(3.33) 


-£>2n 

2 n 


(2 n)\ 
6 n 2 2n 


n 


E <-)' +1 


^ p(n+l-0/ \ 

(/-!)! [) 


(2 n)\ 
6 n 2 2n 


J2n l+1 p {n+l ~ l \n) 

i=i 


( n ) 

a) . 


The coefficients for 1>2 become 



(2 n)! 

6 n 2 2n 


E<-)' +2 

1=2 


p(n+l-l)( 




a 


(n) 


.Mi Vr-V+" ptn+r-o^ 

6 n 2 2n ' (/ _ 1)! 1 

(2n)! s(n,n) (1) (2n)! 

6 n 2 2n ( n _ i)f 1 J 6 n 2 2ri n! ’ 


(3.34) 


(3.35) 


where s(l,k) denote the Stirling numbers of first kind l('32l I33l, [34ll . Formula (3.32) can be com¬ 
puted for different k-values. As an example we present here the computation for k=0. The 
corresponding calculation for k=l and 2 is shown in the appendix section. We find for k=0: 


Example 3.5 


n r)2nan 

B-)' +I (3.36) 

or 

n o2 nan 

E(-) ,+I pV, w = • <3J7) 

with 

n n 

J2(~) l+1 P in+1 ~ l) (n) = (-) n+1 J2(~) l+1 p(l \ n ) ■ (3.38) 

i=i i=i 

4 Application to a series representation of (( 3 ) 

As mentioned in the introduction our polynomial representation of the Bernoulli numbers can be 
used to compute rather fast converging sequences, for example, for C(3). Using (1.1) and (3.39) 
it follows 


1 

C(2 n) = -(2n)C(2)”£ (-)' +1 P“>(n). 

1=1 


(4.1) 
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(4.2) 


In other words this special P-polynomial sum converges rather fast as n increases 


lim (—) l+1 P^(n) = 0 . 

71—»• OO ' ^ 

1=1 


This is essential because the prefactor in (3.39) growths much faster with n as the Bernoulli 
numbers do. As a consequence the P-polynomial sum must compensate this behavior so that the 
product results in B 2n . Also, a representation of C(2n) in terms of £(2) is possible with the help of 
the polynomial representation. Furthermore, by use of (5.1) a direct representation of the In sin(x) 
function is possible. We find 


In sin(irx) = ln(irx) — 2^^ —— 

71= 1 

oo n 

= ln{Trx) - J^C(2) n ^ (-) l+1 P {l] (n) x 


71= 1 


1 = 1 


Following Bendersky [26] we have 

C(3) = -6tt 2 


7T 


'0 Jo 


In sin(ny)dy — — ln( 2) . 


By integrating (5.3) twice and combining the result with (5.4) it follows 


= T “ u ln (z) + tE 


2 oo 


(-U 


2 n 


3 ^2n(2n + 1)(2??, + 2) \6v^6 J 

__ x 2 n 

20 E (2n + 1 )(2n + 2) V6V67 


7T 


^ ^ + i? 


7T 


(4.3) 


(4.4) 


6^6 J 


1 .(4.5) 


If we consider the first two sums only £(3) follows with an error 5 < 10 - '. The convergence is 
very fast, as higher order sums according to the polynomials p6)(n) with index i = 3,4,... start 
with a summation index n = i instead of n = 1. A further advantage of this type of summation is 
that all sums appearing can be expressed in terms of elementary functions based on logarithmic 
expressions. This allows a more detailed insight on £(3) and in consequence on all other zeta- 
values evaluated at odd integer numbers. In the next section we present an alternative series for 
C(3), C(5) and C(7) which converge even faster using Eq. (3.33) for k=2. 


5 Computation of ((3), ((5) and £(7) 

As the main application we use now an alternative polynomial representation to compute £(3), 
C(5) and ((7). With Eq. (3.33) it follows: 

Proposition 4.1 


In sin(nx) 


ln{nx) 

2((2) n //_y+i fn + 2- 

ik ^~ 1)2 n{^ { 2 



(5.1) 
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As a consequence we find for £(3) by integration 


with 


CP) = £ -tlnQ +36 f:(-r^f^) : , 


12 V 3 


n= 1 


2 n 


6^/6 J 


Cl = 


E 


n(n + l)P^(n) f i r \ 2n 


—' (2 n - l)2n(2n + 1)(2n + 2) ^6-^/ 




n(n + l)P (2 )(n + 1) / 


E /H HT ' 'l»T l| / 7T 

n=i (2n + l)(2n + 2)(2n + 3)(2n + 4) V6^6 


2n 


* = E 


n(n + l)P( 3 \n + 2) / 


7T 


—^ (2n + 3)(2n + 4)(2 tt, + 5)(2tt. + 6) \6v^6 


n=1 


2 n 


c 4 = 


E 


n(n + l)P( 4 )( n + 3) / 


7r 


^ (2rz + 5)(2n + 6)(2 tz + 7)(2tt. + 8) \6\/6 


2 n 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 


and this way for higher coefficients c n . Using these first four coefficients only we find the follow¬ 
ing approximate value for £(3): 


C(3) > 1.2020569031595738..., 


(5.7) 


with an error 5 < 0.2* 10~ 13 . The convergence is about 3 orders of magnitude with the summation 
index n. This is about 1 order of magnitude faster compared to the corresponding numerical values 
which result from the series representation (1.7). Furthermore, this new series shows up with an 
alternating sign which has some benefit in estimating its the convergence properties. A further 
advantage of this type of summation is that all infinite sums appearing can be expressed in terms 
of elementary functions based on logarithmic expressions. This allows a more detailed insight on 
C(3) and in consequence on all other zeta-values evaluated at odd integer numbers. The explicit 
comparison is shown in Table III: where the numerical comparison has been performed between 
our approach and the BBP-formalism by use of the following formula If35l : 


57r 2 /n(3) 3ln(3Y 


1 


104 


104 


E 


ALY 


729 


+ 


243 


81 


1053 f^V 729 / V(12n + 1) 3 (12n + 2) 3 (12n + 4) 3 

54 27 


n =1 
81 


(12n + 5) 3 (12n + 6) 3 (12ra + 7) 3 

9 3 3 


(12n + 8) 3 (12n + 10) 3 (12n + ll) 3 (12ra + 12) 3 

Analogously we find for £(5): 


(5.8) 


3vr 2 

C(5) = jjCP)- 


7r 


12528 1044 Y3- 


1447T 2 
29 


£(-r + v 


( 7T \ 


2 n 




n= 1 


\6VqJ 


(5.9) 
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C(3) 

C(3)-Zeta series 

C(3)-(Zeta series+BP) 

C(3)-BBP formula If35l 

1st order 

n=l 

<5=0.2* lO" 04 

5=0.2* 10~ 05 

5=0.5* 10 -06 

2nd order 

n=2 

5=0.2* lO" 06 

5=0.1*10~° 8 

5=0.9* 10" 10 

3rd order 

n=3 

<5=0.3* 10" 08 

5=0.8* 10" 12 

5=0.4* lO' 13 

4th order 

n=4 

<5=0.4* 10~ 10 

5=0.2* 10" 13 

5=0.2* lO' 16 


Table 3: Approximate computation of £(3) as a function of the summation index n by use of 
Eq. (1.7) without and with use of the polynomial representation. The numerical errors are com¬ 
pared to the BBP-type formula Il35l 


with 


Cl = 


E 


n(n + l)-P^(n) 


^ (2 n - l)2n(2n + l)(2n + 2)(2n + 3)(2n + 4) \6V6 J ’ 


A) 


2 n 


(5.10) 


*= E 


n 


(n + l)P( 2 \n + 1) 


(2 n + l)(2n + 2)(2n + 3)(2n + 4)(2n + 5)(2n + 6) \6VQ J 


c 3 = 


n=1 


E 


A) 


2 n 


n(n + 1 )P( 3 \n + 2) 


^ (2ri + 3)(2n + 4)(2n + 5)(2n + 6)(2n + 7)(2n + 8) ^6^6 


7r 


2 n 


(5.11) 


(5.12) 


c 4 = 


E 


n(n + l)P (4 ^(n + 3) 


(-Cl 


2 n 


(2 n + 5)(2 n + 6)(2ra + 7)(2ra + 8)(2ra + 9)(2 n + 10) \6V6J ’ 


(5.13) 


Table IV shows the corresponding computational results: where the numerical comparison has 
been performed between our approach and the BBP-formalism by use of the following formula 

( 3611371 : 


31 

32 


C(5) 


343 4l , . 5 

-7r 4 /n(2) H- 

99360 v ; 2484 

OO 


7r 2 (/n,(2)) 3 - —(Zn(2)) 


128 

"69" 


E 

n= 1 


J_\ n cos(!f) 

V2 J 


nr 


_ 20 y. 


n= 1 


1035 
1 
2 


rr 


(5.14) 
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C(5) 

C(5)-Zeta series 

C(5)-(Zeta series+BP) 

C(5)-BBP formula EH 

1st order 

n=l 

5=0.6* 10 -06 

5=0.1* 10~ 06 

5=0.7* lO" 06 

2nd order 

n=2 

5=0.4* 10" 08 

5=0.1*10~° 9 

5=0.1*10 -08 

3rd order 

n=3 

5=0.3* 10- 10 

5=0.1* 10“ 12 

5=0.8*10“ 11 

4th order 

n=4 

5=0.3* 10- 12 

5=0.1* 10“ 15 

5=0.1*10“ 12 


Table 4: Approximate computation of £(5) as a function of the summation index n by use of 
Eq. (1.15) without and with use of the polynomial representation. The numerical errors are com¬ 
pared to the BBP-type formula 071 


For £(7) we found: 


72t r 2 2vr 4 


49vr 6 


7T , /7T 


^ 659 ^ 1977^ 3 ' ) + 5337900 + 266895““ V3 


In 


34567r 4 g. ( )w+1; 


5931 


72 — 1 


A) ■ 


with 


Cl = 


E 


n(n + l)pW(n) 


272 


c 2 = 


E 


(2n - l)2n(2n + l)(2n + 2)(2n + 3)(2n + 4)(2n + 5)(2n + 6) 

/ \ 272 

n(n + l)p( 2 )(n + 1) (^J 


(5.15) 


(5.16) 


(2n + l)(2n + 2)(2n + 3)(2n + 4)(2n + 5)(2n + 6)(2n + 7)(2n + 8) 


(5.17) 


c 3 = 


E 


n(n + 1 )P( 3 \n + 2) 


6^6 


272 


(2n + 3)(2n + 4)(2n + 5)(2n + 6)(2n + 7)(2n + 8)(2n + 9)(2n + 10) 


(5.18) 


OO 


n(n + l)p( 4 )(n + 3) 


6^ 


272 


C4 ( 2 n + 5) (2 n + 6)(2n + 7)(2n + 8)(2ra + 9)(2ra + 10) (2 n + 11) (2n + 12) 


(5.19) 


Table V shows the corresponding computational results: where no BBP-type formula is available 
for C(7). 
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C(V) 

C(7)-Zeta series 

C(7)-(Zeta series+BP) 

C(7)-BBP formula 

1 st order 

n=l 

5=0.7* KT 08 

5=0.2* 10“ 08 


2nd order 

n=2 

5=0.3* lO" 10 

5=0.1*10~ n 


3rd order 

n=3 

5=0.2* 10- 12 

5=0.8* 10" 15 


4th order 

n=4 

5=0.2* IQ' 14 

5=0.6* 10" 18 



Table 5: Approximate computation of £(7) as a function of the summation index n by use of 
Eq. (1.16) without and with use of the polynomial representation. For ((7) no BBP-type formula 
is available for reasons of comparison. 


6 Approximations of zeta numbers using B 2n and B 2n _2 

We present a formula which allows to compute the Bernoulli number B 2n as a function of the 
Bernoulli number B 2n -2 only. This formula decouples B 2n from all other Bernoulli numbers. We 
simply compute (3.32) for k=4 with the help of (3.3). It follows: 

-i „ 


B- 


2 n 


2 n 


B 2n - 2 
2n — 2 


+ (-) 


ra+l 


(2 n)\ 
2 2 n 6 r 


2 n 
4 




l+l 


1=1 


n + 4 — l 
4 


P {l) ( 


n ) 


( 6 . 1 ) 


This formula allows the computation of B 2n as a function of B 2n _2 only by use of the correspond¬ 
ing Bernoulli-type polynomials P+n). Furthermore, it allows for an approximation of the odd 
zeta numbers which works faster by about one order of magnitude when compared to (5.2). It 
follows first: 


71+1 




£+1 


With this we can write: 


i=i 


n + 5 — l 
4 


P {l \n + 1) - 


2n(2n + 1) 
4-7T 2 


C(2n + 2) . 


C(3) = 


7r 


+ 


7T 


15552 8 


- —In 
12 


7r 


+ 


7T 


2 oo 




71+1 


(. 


7r 


n= 1 


with 


Cl = 


\6V6 

n + l)(n + 2)(n + 3)(n + 4)P^)(n + 1) 

24 L ^ (2n - l)2n(2n + 1)(2 n + 2)(2 n + 3)(2n + 4) 

1 3 


( 6 . 2 ) 


(6.3) 


E 


24n(n + l)pW(n) 


(2n — l)2n(2n + l)(2n + 2)(2n + 3)(2n + 4)J \Q^/QJ 


2n 


(6.4) 
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C 2 = 


E 


n(n + l)(n + 2 )(n + 3 )P^ 2 \n + 1) 


24 L^ (2 n- l)2n(2n + l)(2n + 2)(2ra + 3)(2ra + 4) 

24„(„+l)pW(»+l)( s 5-) 2 Ujl_\ 


(2 n + l)(2n + 2)(2n + 3)(2n + 4)(2n + 5)(2n + 6)J \6V6J 


2 n 


(6.5) 


c 3 = 


_■ OO 

a[E 


n(n + l)(n + 2 )(n + 3 )P®(n + 2) 

^ (2n + l)(2n + 2)(2n + 3)(2n + 4)(2n + 5)(2n + 6) 

24n(n +l)P (3) (n + 2) 

(2rz + 3)(2n + 4)(2ra + 5)(2ra + 6)(2 n + 7)(2n + 8)J \6v^ 


2 n 


( 6 . 6 ) 


c 4 = 


OO 

[E 

n =1 


n(n + l)(n + 2)(n + 3)P^(n + 3) 

24 L ^ (2ra + 3)(2n + 4)(2n + 5)(2n + 6) (2 n + 7) (2 n + 8) 

24 n ( n +l)F<‘)( n + 3)( 5 5 s ) 2 ,/^N 

(271 + 5)(2ra + 6)(2ra + 7)(2ra + 8)(2ra + 9)(2n + 10)J / 


2 n 


(6.7) 


and this way for higher coefficients c n . Using these first four coefficients only we find an ap¬ 
proximate value for £(3) with an error 5 < 0.4* 10 -14 . The convergence is about 4 orders of 
magnitude with the summation index n. This procedure allows for a systematic improvement of 
the convergence behavior if higher values k=6 or k=8 are used in (3.2). As an example we present 
the corresponding formula for k=6: 


2n \ / Ehn _ 5 B 2n -2 1 P 2 n—4 \ , w + i (2 n)\ \r^,_y + i f n + 6 - 

6 )\2n 2n-2 2n-AJ y ’ 2 2n Q n 1 l 6 


P w (n) . 


( 6 . 8 ) 


Our results establish a new and very fast option to compute zeta and related numbers. As an 
outlook, it could be of great interest to combine our formalism with the BBP approach by use of 
corresponding polylogarithmic identities to further improve the convergence properties in explicit 
computations. 


7 SUMMARY 

In summary, we have presented a unique computational scheme for the explicit calculation of the 
Riemann £ function and its first derivatives at all positive and negative integer values. This way 
we have shown that all these numbers are directly attributed to Bernoulli numbers, but with an 
increasing level of complexity when going, for example, from £(2) t0 a related sub-sum like U 4 . 
The computational scheme is based on a new polynomial representation of the Bernoulli numbers 
in connection with Bendersky’s L-numbers, which appear in context with the logarithmic Gamma 
function. As a first application we performed approximate calculations of £(3), C(5) and Q(l) in 
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terms our polynomial representation, where this computational procedure is applicable to all £- 
values with integer arguments, as well as to related numbers like Catalan’s constant. Finally, 
we have shown that a computation of B 2n as a function of B 2n -2 or, for example, of B 2n -2 and 
B 2 n ,—4 only is possible by use of the polynomial representation. The result is a further improved 
approximate computation of these numbers. 


8 APPENDIX A 


In analogy to example (3.5) for k=0 it results for k=l and k=2: 

n ( o2n 0\R n 

J2(~) l+11 P (n+1 ~ l \ri) = ° 

i=i 


B- 


2 n j 


or 


with 


^(-) z+1 /P w ( 


n = — 


\n+l 


1=1 


2n(2n\) 

(4 n -{n- l)2 2n )6 n 
2n(2n!) 


B- 


2 n 


^(-) /+1 / p (n+1 


l ^n) = 


i=i 


, n+1 2n(2 2 " — 2) 


4n 


(n - 1)2 2 " ^ 


£(-)' +1 l P ®( n ) . 


For k=2 we find 


or 


B-r ? n) = P+(2 ? -W ^ [ 


«=1 


(2d) 


^(-)* +1 Z 2 P (z) (n) 


2=1 


(-) 


n+1 (4n(2n + 3) + (n 2 - 6n + l)2 2n + 


(2n!)2n 


B- 


2nt 


with 


J2(~) l+1 l 2 P {n ~ l+1 \ri) = 


i=i 


(-) 


n+1 


2n(2 + (2n - 2)2 2n ) 


4n(2n + 3) + (n 2 — 6 n + l)2 2n 

where the simple symmetry property 


X>)‘ +1 i 2 P { '\n) , 


^2(-) l+1 r p w ( 


n = 


i=i 


i=i 


( 8 . 1 ) 


( 8 . 2 ) 


(8.3) 


(8.4) 


(8.5) 


( 8 . 6 ) 


(-) n+1 ^(-)* +1 (n + 1 - l) m P in+1 ~ l \n) , (8.7) 


which we have used before in the case of m=l has been applied here again for m=0,l and 2. 
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9 APPENDIX B 


Here we show the Bernoulli polynomials up to n=6. It follows: 


P (5) (n) 


P^iyi) = — 
n 

(9.1) 

p,2) <"> = 2 ! 5 

(9.2) 

P m (n) _ 3< 21 ™ - «) 

^ 2 3 * 5 2 * 7 

(9.3) 

. 63n 2 — 387n + 590 

P W= 2 4 * 5 3 * 7 

(9.4) 

3(4851n 3 - 59598n 2 + 242737n + 327210) 

2 7 5 * 4 *7 2 * 11 

(9.5) 


P (6) (n) 


3(189189n 4 - 3873870n 3 + 29616015n 2 - 100104550n + 126087736) 

2 8 * 5 6 * 7 2 * 11 * 13 


(9.6) 
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